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Analysis of N* spectra using matrices of correlation functions based on 
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We present results for ground and excited-state nucleon masses in quenched lattice QCD using anisotropic 
lattices. Group theoretical constructions of local and nonlocal straight-link irreducible operators are used to obtain 
suitable sources and sinks. Matrices of correlation functions are diagonalized to determine the eigenvectors. Both 
chi-square fitting and Bayesian inference with an entropic prior are used to extract masses from the correlation 
functions in a given channel. We observe clear separation of the excited state masses from the ground state mass. 
States of spin > 5/2 have been isolated by use of G2 operators. 



1. INTRODUCTION 

Reproducing the spectrum of baryon reso- 
nances with spin-1/2 and spin-3/2 and both pari- 
ties is an important test of lattice QCD. For that 
we require three-quark operators that transform 
irreducibly under the spinorial rotation group of 
the lattice 0. Local and nonlocal straight-link 
operators corresponding to irreducible represen- 
tations (IRs) Gig^u, Hgu and G2g.u are con- 
structed to obtain suitable sources and sinks [5]. 
Here we analyze N* spectra using these opera- 
tors. 

To determine the N* excited states, a matrix of 
correlation functions is computed in the quenched 
approximation to QCD using irreducible baryon 
interpolating operators Bi{x,t) of definite quan- 
tum numbers, 



The Cij{t) matrices for Gi and H states for 
each parity are constructed using local operators, 
smeared local operators and smeared straight-link 
operators. For G2, we have one operator which is 
of smeared straight-link type (Tables 1, 3 in 

The computation of masses of the lowest-lying 
resonances is based on the variational method ap- 
plied to the matrix of correlation functions. In 
this paper we solve the generalized eigenvalue 
equation, 

C^J (t) f/") (t) - A(") {t, to) a, (to) f/") (t) (2) 

and determine eigenvectors I^'^"^(i) for each t, 
with to close to the source time. Then the masses 
of N* states correspond to the eigenvalues of 
Eq. ©: A(")(t,io)-t»t;? e-™'°'(*-*o) 3. The ef- 
fective masses are determined from 
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Another way to extract spectrum information 
is to calculate the spectral mass density p{uj) from 
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Figure 1. G\g effective masses for a selected few 
low-lying states. 
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Figure 2. Lowest positive parity effective masses 
for the IRs Gig, Hg and G2g- 



lattice correlation functions using the Maximum 
Entropy Method (MEM) g], 

C{tM) — y dc.p(c.)e^"(*-*°). (4) 

One of the advantages of MEM is that it utilizes 
data on a wide range of available time slices and 
has been shown to yield results even for noisy 
data g]. This feature may be helpful in extracting 
masses of excited states. 

2. RESULTS 

We use an ensemble of 287 quenched, 
anisotropic 16^ x 64 lattices with renormalized 
anisotropy ^ = 3.0 and /3 = 6.1, corresponding to 
a^^ = 6.0 GeV We use the anisotropic Wil- 
son action. The parameters of the Wilson fermion 
action are tuned nonperturbatively so as to sat- 
isfy the continuum dispersion relation E{p)'^ = 
E{0)'^ + c(p)2p2 at a pion mass ~ 500 MeV. 
To improve the coupling of operators to the lower 
mass states we employ gauge-covariant smear- 
ing of the quark fields on both source and sink: 
ijj{x) = (1 + a^A{U)/AN)'^iP{x), where A(U) is 
the three dimensional Laplacian and U denotes 
APE-smeared SU{3) link variables. The param- 
eters used to smear the quark fields are cr = 3.6 
and TV = 32. 

The effective masses are calculated from 10 x 10 



Gig matrices with to = 2. In Fig. we choose to 
show a few low-lying states that are clearly sep- 
arated. However, the details of the states above 
the ground state are under study. The plot shows 
a good plateau for the ground state and statisti- 
cally significant splittings for a couple of excited 
states. 

In Figs. 12 and 13 we have collected the effective 
mass plots of the lowest states of both parities 
for Gi, G2 and H. The ratio of lowest masses 
for Gig and Gi„ is roughly in accordance with 
experiment, for spin-1/2 states, the Gi„ mass be- 
ing higher. The effective masses for Hg^u are ob- 
tained using a 7 X 7 matrices of correlation func- 
tions. The lowest negative parity Hu state has 
smaller mass than the lowest positive parity Hg 
state. This is compatible with the pattern found 
in nature for spin-3/2. However, the masses of 
Hu and Gi„ overlap within errors. Our result for 
G2 masses also reveals reasonable separations of 
the G2g and G2U masses, G2g being lower. From 
Fig. 13 it is evident that the effective mass for 
Hg (allowed spin 3/2+, 5/2+, ■ • •) is very similar 
to that for G2g (allowed spin 5/2+, 7/2+, • ■ •)• 
These states are orthogonal. One possibility for 
this is that the lowest Hg state has spin-3/2+ and 
its mass is accidentally close to that of the lowest 
G2g state. Another possibility is that the low- 
est Hg state is spin-5/2+, in which case the same 
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Figure 3. Lowest negative parity effective masses 
for the IRs Gi„, iJ„ and G2u- 

state must be present in Hg and G2g, but not in 
Gig. Study over different values of lattice spacing 
is required to decide. 

Finally, we present the G23 MEM spectral func- 
tion in Figure 01 We find that the peak of the 
spectral density roughly corresponds to the effec- 
tive mass value. 

In Tabled we summarize our preliminary esti- 
mates of the lowest masses for the different repre- 
sentations extracted from single-exponential fits 
to the X''°'\t,to) of Eqn.[l] The effective masses 
for the lowest states of Gi, H and G2 for both 
parities, whether obtained from the variational 
method or preliminary MEM analysis, show a 



Table 1 

Numerical estimates from single-exponential fit to 
A''"-'(i,to) corresponding to the lowest masses of 
N* spectra. 



IRs 


Fit range 


"loft 


~ Mass (MeV) 


Gig 


9-20 


0.208 (4) 


1250 


Glu 


8-12 


0.321 (4) 


1930 


Ha 


9-12 


0.410 (2) 


2460 


Hu 


6-15 


0.315 (4) 


1890 


G2g 


9-14 


0.409 (7) 


2450 


G2U 


8-15 


0.475 (7) 


2850 




Figure 4. Example of a spectral density function 
from a MEM analysis of G2g correlation func- 
tions. 

spectrum of distinct N* masses. However, the 
behavior of the spectrum with 711^^ and the sensi- 
tivity of the spectrum to variations in the lattice 
volume has yet to be studied. 

This work is supported by US National Sci- 
ence Foundation under Awards PIIY-0099450 
and PHY-0300065, and by US Department of 
Energy under contract DE-AC05-84ER40150 and 
DE-FG02-93ER-40762. 

REFERENCES 

1. C. Morningstar et ai, Nucl. Phys. B (Proc. 
Suppl.) 129 (2004) 236 and these proceedings 

2. I. Sato et al. Nucl. Phys. B (Proc. Suppl.) 129 
(2004) 209 and these proceedings 

3. C. Michael, Nucl. Phys. B259 (1985) 58; M. 
Liischer et ai, Nucl. Phys. B339 (1990) 222; 
C.R. Alton et ai, Phys. Rev. D47 (1993) 5128 

4. H.R. Fiebig, Phys. Rev. D65 (2002) 094512 

5. R.G. Edwards et ai, Nucl. Phys. B (Proc. 
Suppl.) 119 (2003) 305 



